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Abstract 

In this paper, we study the asymptotics of the discrete Chebyshev polynomials t n (z,N) as the 
degree grows to infinity. Global asymptotic formulas are obtained as n — > oo, when the ratio of the 
parameters n/N = c is a constant in the interval (0, 1). Our method is based on a modified version of 
the Riemann-Hilbert approach first introduced by Deift and Zhou. 



1 Introduction 



< 

u 

The discrete Chebyshev polynomials were introduced by Chebyshev in the study of least-squares data 
fitting. They are explicitly given by 

> . t n (z,N)=n\A n i [ Q( Z - n N )}, n = 0,l,..-,N-l, (1.1) 
vo " 

£^ | where A is the forward difference operator with unit spacing on the variable z; see [14, (2.8.1)]. These 

{SJ " polynomials are orthogonal on the discrete set {0, l,--- ,N — 1} with respect to the weight function 

■ w(k) = 1, k = 0, 1 • • • N — 1. More precisely, we have 
O" 

^ ^ (N + nV 

Uk, N)t m (k, N) = — - ^ _ } n _ — 6 n , m , n,m = 0,1,- •• ,N-1. (1.2) 

Like all other classical orthogonal polynomials, these polynomials also satisfy a three-term recurrence 
" relation. Indeed, we have 

(n + l)t n+1 (z, N) = 2(2n + l)(z - ±(iV - l))t n (z, N) - n(N 2 - n 2 )t n _!(z, N); (1.3) 

see Gautschi [6]. 

The discrete Chebyshev polynomials are also known as a special case of the Hahn polynomials [2, 
p. 174] 

Q n (z; a, (3, N) := zF 2 {-n, -z, n + a + + 1; —N, a + 1; 1); (1.4) 
see also [8]. Recall the weight function of the Hahn polynomials 

w(z;a,P,N) := — — - _ ^ , z = 0, 1, • • • , N - 1. (1.5) 
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Note that by taking a = ft = 0, the weight function becomes 1. In this case, the Hahn polynomials reduce 
to the discrete Chebyshev polynomials except for a constant factor. 

For further properties and applications, we refer to Szego [14] and Hildebrand [7]. Like other dis- 
crete orthogonal polynomials, these polynomials do not satisfy a second-order linear differential equation. 
Hence, the asymptotic theory for differential equations can not be applied. By using the steepest descent 
method for double integrals, Pan and Wong [12] have recently derived infinite asymptotic expansions for 
the discrete Chebyshev polynomials when the variable z is real and the ratio n/N lies in the interval 
(0, 1); their results involve a confluent hypergeometric function and its derivative. For more information 
about the integral methods, we refer to Wong [17]. On the other hand, Baik et al. [1] have studied the 
asymptotics of discrete orthogonal polynomials with respect to a general weight function by using the 
Riemann-Hilbert approach. The results in [1] are very general, but it is difficult to use them to write out 
explicit formulas for specific polynomials. In the case of the Hahn polynomials (1.4), only the situation 
when both parameters a and ft are large has been considered. More precisely, the authors of [1] considered 
the case when a = NA and ft = NB, where A and B are fixed positive numbers, thus excluding the 
special case of the discrete Chebyshev polynomials (i.e., a = ft = 0). Moreover, the results in [1] are more 
local in nature; that is, more formulas are needed to describe the behavior of the orthogonal polynomials 
in the complex plane. 

In this paper, we investigate the asymptotics of the discrete Chebyshev polynomials as n — > oo, when 
the ratio of the parameters n/N is a constant c £ (0, 1). Global asymptotic formulas are obtained in the 
complex z-plane when n goes to infinity. Our approach is based on a modified version of the steepest- 
descent method for Riemann-Hilbert problems introduced by Deift and Zhou in [4]; see [3], [11] and [15]. 
More precisely, we derive two Airy-type asymptotic formulas for the discrete Chebyshev polynomials. The 
regions of validity of these formulas overlap, and together they cover the whole complex plane. 

The presentation of this paper is arranged as follows: In Section 2, we first present the basic Riemann- 
Hilbert problem associated with the discrete Chebyshev polynomials. Then, we transform the basic 
Riemann-Hilbert problem into a continuous one, which is similar to that in [3, 15]. In Section 3, we 
introduce some auxiliary functions, known as the ^-function and the (^-function. Moreover, we define a 
function D(z), which is analogous to the function first used in [15] to remove the jumps of the continuous 
Riemann-Hilbert problem near the endpoints of the interval of orthogonality. In Section 4, we construct 
our parametrix by using Airy functions and their derivatives, and prove that this parametrix is asymp- 
totically equal to the solution of the Riemann-Hilbert problem associated with the discrete Chebyshev 
polynomials formulated in Section 2. Our main result is given in Section 5. In Section 6, we compare our 
formulas with those given by Pan and Wong [12]. 



2 Riemann-Hilbert Problem 

We start with the fundamental Riemann-Hilbert problem (RHP) for the discrete orthogonal polyno- 
mials. Note that the leading coefficient of t n (z,N) is (2n)!/n! 2 . Hence, the monic discrete Chebyshev 
polynomials are given by 

n' 2 
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Clearly, they satisfy the new orthogonality relation 

N-l 

(k)lVN,m(k) — C>n,m./77V,ra; 



(2.2) 



where 



fc=o 



2 (2n + l)r(iV-n)r 2 (2n + l) 

77V,n =~ 



T(iV + n+l)r 4 (n + l) 
Consider the following RHP for a 2 x 2 matrix- value function Y{z): 

( Y a ) Y(z) is analytic for z G C \ N; 

( Yj,) at each z = k S N, the first column of Y is analytic and the second column of Y has a simple pole 
where the residue is 

ResYO) = ]xmY(z) 



z=k 



( Y c ) as z — > oo, 



Y(z) 



[ - 

Z 











z~ n / 

By a well-known theorem of Fokas, Its and Kitaev [5] (see also Baik et al. [1]), we have 
Theorem 2.1. The unique solution to the above RHP is given by 



( 



Y{z) :-- 



7T N>n (k) \ 



k=0 



z — k 



2 N-i niff^mn-i(k) 

TN,n-l^N,n-l{z) J2 ~ 
k=0 



z — k 



(2.3) 



Let X]y denote the set defined by 



X N := {x N)k }^ = Q, where x Njk :-- 



k + 1/2 
N ' 



The »i\rfc' s are called nodes and they all lie in the interval (0, 1). Following Baik et al. [1], we introduce 
the first transformation 



H(z) 



N~ n 
N' 1 



n Y(Nz - 1/2) 



/N-l 

\[{Z- Xnj)" 1 
3=0 



\ 









N-l 

n i z - x N,j) 

3=0 



N- 



Y(Nz- 1/2) 



N-l 

n ( z - x nj) 

3=0 



(2.4) 



where o~ 3 : = 



1 



is a Pauli matrix. 



A straightforward calculation shows that H(z) is a solution of the RHP: 
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H a ) H(z) is analytic for z £ C \ X/v; 



( Hb) at each x^,k G -XjV) the second column of H(z) is analytic and the first column of H(z) has a simple 
pole where the residue is 



Res H(z) = lim H(z) 



( o 

N-l 







n {z-xMj)- 2 o 



3=0 



( H c ) as z — » oo, 



H(z)=(l + 0(- 



y n-N 







-n+iV 



We next consider the second transformation which will remove the poles as well as transform the 
discrete RHP into a continuous one. Let 5 > be a sufficiently small number, and we define 



R(z) := H(z) 



1 



\ 



for z G O-i- , and 



Tie ±Nmz cos(Nirz) 



R(z) := H(z) 



(2.5) 



(2.6) 



for z ^ Q±, where the domains Q± are shown in Figure 1. Let E_|_ be the boundary of f2+ in the upper 
half-plane, and £_ be the mirror image of E+ in the lower half-plane. For the contour E = (0, 1) U E±, 
see also Figure 1. 




1 + iS 



1-iS 



Figure 1: The domains Q± and the contour E. 



Lemma 2.2. For each & S X^v, the singularity of R(z) at x^ k is removable; that is, Res R{z) = 0. 

Z=%N,k 
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Proof. For z £ tt±, it follows from (2.5) that 



Rn(z) = H u (z) + H 12 (z) J ^ y R 12 (z) = H 12 (z). (2.7) 



Since H\ 2 {z) is analytic by (Hf,), the residue of R\ 2 (z) at Xjvft is zero. From we observe that the 

residue of Hu{z) at z = x^,k is 



N-l 

Res iJn(^) = Hi 2 (x N k ) TT (xAr fc - x N j )~ 2 . (2.8) 

2=Z]V,fc 

On the other hand, it is readily seen that 



3=0 



N-l 



Res J 1 L = _ JJ ^ XJV . -2_ 2>9 

Thus, applying (2.8) and (2.9) to (2.7) shows that the residue of Ru(z) at z = xjy^ is zero. The entries in 
the second row of the matrix R(z) can be studied similarly. This completes the proof of the lemma. □ 

Note that this transformation makes R+(x) and R-(x) continuous on the interval (0,1). Therefore, 
R(z) is a solution of the following RHP: 

(R a ) R(z) is analytic for z G C \ E; 

(Rb) the jump conditions on the contour E: for x £ (0, 1), 

R+(x) = R„(x) ( f), (2.10) 

where 



for z E E±, 



where 



r(x) 1 J 



, , 4cos 2 (iV7ra;) 
r ^ = N-l ' ( 2 - U ) 

i=o 



R + (z)=R.(z){^)^ °), (2.12) 



(i? c ) as z — > oo, 



r±(z) 1J> 

r ±(z) = T^p& ; (2.13) 

^ n (z-*n,j) 2 

j=0 



i-))( Zn N z -n + N 



For simplicity in the following section, here we have introduced the notation W := 2N?ri. 
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3 The Auxiliary Functions 

For our subsequent analysis, we need some auxiliary functions. Consider the monic orthogonal poly- 
nomials PN,n(z) ■= N~ n -KN tn (Nz — 1/2). From (1.3) and (2.1), we have 

PN,n+l(z) = (z- a n )p N:n (z) - 6 n p7V jn _i (z) , (3.1) 

where 

1 1 - (n/N) 2 

a n — 77; b n — —— — -. 

2 16 — A/n z 

Using the method introduced by Kuijlaars and Van Assche [9], we can derive the equilibrium measure 
corresponding to the discrete Chebyshev polynomials in our case. Recall c := n/N, and let 

a:=\-\^J\-c 2 , + (3.2) 

The density function is given by 



— arcsinf ° - ), xe[a,b], 
vrc 2yx - x 2 

-, x £ [0,a] U [b, 1]. 
^ c 



(3.3) 



As usual, we now define the auxiliary g-function and ^-function. 
Definition 3.1. The g-function is the logarithmic potential defined by 



g(z):= f \og{z-s)»(s)ds, zeC\(-oo,l], (3.4) 

Jo 



and the so-called 4>-function is defined by 



m ■■= \ - 9{z) (3.5) 



for z € C \ (—00, 1], where I := 2 log \b — s\fi(s)ds is called the Lagrange multiplier. 

An explicit formula for the g-function is given in (6.25) in Section 6. On account of (3.3), the derivative 
of g(z) is given by 

n s 2, / z+ J(z-a)(z-b) + c/2 \ 1 fz-l\ , n e , 

g'( Z ) = - log V V A J 1 + - log . 3.6) 

c \(z-l) + y/(z-a)(z-b)-c/2j c & \ z J 

Introduce the ^-function 

0*( z ):=^ (-g'(s) T ^irij ds = <f>(z) ± - z) (3.7) 
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for z G C±. It is readily seen that 



W A U S I 1-8- y/(s-a)(s-b) + c/2 J c & \ s 



for z G C\(-oo,6] U [l,oo). 
Similarly, we also set 



[-g'(s)T^mjds = (j)±m(l-^z) (3.9) 



for z G C±. Note that for x G (0, a), we have from (3.4) 

9±( x ) = / log \x — s|/i(s)<is i 7ri / fi(s)ds. 
Jo Jx 

Thus, from (3.9) and (3.5), it follows that </>+(x) = <f>-(x) for x G (0, a); that is, <p(z) can be analytically 
continued to the interval (0, a). 

The functions <f>(z) and (f>*(z) play an important role in our argument, and the following are some of 
their properties. 

Proposition 3.2. For x G (a,b), we have 

<t>*+(x) + <l>*_(x) = 0, 4> + {x) + 4>-(x) = 0. (3.10) 

Furthermore, we have 

0*(x)<O /orxe(M) « nd #c)<0 /orxG(0,o). (3.11) 
For any x G (6, 1) and sufficiently small 5 > 0, we have 

Recj)*(x±i5)=4>*(x) + 0(5 2 ), Rc<f>(x ±i5) = 4>*(x) - -tt5 + 0(5 2 ). (3.12) 

c 

For any x G (0, a) and sufficiently small 5 > 0, we have 



Recf(x±i6) = cf){x) + 0(5 2 ), Rec/>(x ± iS) = <f>(x) - -ir6 + 0(5 2 ). (3.13) 



Proof. Since 



jl.^ T/ 2 ! ( s±iy/(s-a)(b-s) + c 2 \ 1, A-s\\ , ,„... 
</4 x = - / - log vv ===== ; + - log ds 3.14 

for x G (a, 6), it is easy to verify that 

ft(x) + * W = - jT (flog ( (1 .',; ( r +e) ) + I log (^)) «. = 0. (3.15) 

Similarly, we also have <j>+(x) + <fi-(x) = for x G (a,b), thus proving (3.10). Together with (3.5) and 
(3.7), this implies that log |x — s\/j,(s)ds = 1/2 for x G (a, 6). 



For any small e > and z E U(b, e) := {z E C : \z — b\ < e}, we have from (3.8) 

-O(e) 



n*) = 



2 / V(fe-a)( 8 -ft) , ^/(b-g)( s -&) 



+ 



3c 2 



b + c/2 l-b + c/2 

(l_ c 2)l/4 (z _ 6) 3/2 +0 ( £ 2)_ 



r/.s 



Here, we have used the fact that a + b = 1, 6— a = (1 — c 2 ) 1 / 2 and ab = c 2 /4; see (3.2). By strai; 
calculation, it is readily seen from (3.8) that 



(3.16) 
ditforward 



4>* (x) = — log 
c 



■log 



+ ^{x- a)(x -b) + c/2 



1 - x - y/(x - a)(x - b) + c/2 J c 

( x + y/Jx -a)(x-b) + c/2 
I 1 - x - yj{x - a)(x - b) + c/2 



1 , /l - X 
-log 



1 - X 



< 



(3.17) 



for b < x < 1. Thus, we obtain (f)*{x) < 4>{b + e) < 0. 

In the same manner, if z E U(a, e) := {z E C : \z — a\ < e}, then 



/' 


[!( 


/a 





a + c/2 



+ 



1- a + c/2 



+ 0(e) 



ds 



-^(l-c^l\a-z)^ + 0{^) 



(3.18) 



and 



0'(x) = 



2, / x- V (a- a) (b- s ) + c/2 \ 1, /1-x 

--log — log 

c - x+ s/{a - x){b- x) + c/2 J c 



1 

--log 

c 



+ y/[x - a)(x - b) + c/2 



1 -x 



> 



1 - x - y/(x - a)(x - b) + c/2 J 

for < x < a. Consequently, <f)(x) < <p(a — e) < for < x < a, thus proving (3.11). 
For any x E (6, 1) and sufficiently small 5 > 0, we have from the Taylor expansion 

Recj)*{x ± itf) = Re</)*(x) T o"lm0*'(x) + O(o 2 ) = Re<f (x) + 0(5 2 ). 

The last equality follows from the fact that (j>*'(x) is real; see (3.17). Coupling this with (3.7) 

Re0(x±t<5) = Re^*(x) - -ttJ + 0(<5 2 ); 

c 

thus (3.12) holds. 

Similarly, for any x E (0, a) and sufficiently small 5 > 0, we have 

Re^(s±i£) = Re<Kx) + <nm0 / (x) + O(o 2 ) = Re <£(x) + 0{5 2 ). 

From (3.9), it follows that 



(3.19) 



(3.20) 



gives 



Re^(x±io) = Re0(x) - -tt5 + 0{5 2 ). 



This ends the proof of (3.13). 



(3.21) 

(3.22) 

(3.23) 
□ 
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(a) z-plane 



(b) 0*-plane 




(c) 0-plane 



V . 






<M1) 








Mb) 


u 




M a ) 




Mo) 







(d) 0-plane 



Figure 2: The upper half-plane under the transformations <j>, (j> and <f>* . 

From the definition in (3.9), it is easy to see that 4>(a) = and (f>+(b) = (1 — -)m. Furthermore, from 
(3.8) it follows that 4>*(b) = 0, and from (3.7) and (3.9) we have 4>+( a ) = (~ — l)iri. Let us now give a 
brief outline of the argument used in establishing the following mapping properties of the ^-function. 

Proposition 3.3. The images of the upper half of the z-plane under the mappings (f>* , (ft and <f> are 
depicted in Figure 2. 

Proof. For x € (6, 1), <j)*{x) is real and negative in view of (3.11). Furthermore, (3.17) implies that cp*(x) 
is a monotonically decreasing function in (b, 1). Since (p*(b) = 0, 4>*(1) is negative. For x G (a, b), we have 
from (3.7), (3.5) and (3.4) 

<t>*+(x) = ^ — J log | a; - s\/j,(s)ds + m J ^~ ~ m( s )^ ds. 
On account of a statement following (3.15), the first two terms on the right-hand side cancel. Hence, we 



9 



obtain 



fi(s) I ds 



for x G (a, 6). By (3.3), \ > fi(s) for s G (0,1). Thus, Rec^(x) = and Imc^(x) > for x G (a, 6). 
Furthermore, since (Im</>+(x))' = — 7r(^ — /J,(x)) < 0, Ini(/)^_(x) is monotonically decreasing in (a,b). By 
a similar argument, it can be shown that for x G (l,+oo), Im(/>l(x) is negative and decreasing; also 
Re 0?j_(x) < Re^>?j_(l) and Re<f)+(x) is decreasing. Again, by a similar argument, it can be shown that for 
x G (0,a), Re</>* + {0) =Re^(l), Im^(O) = Im</>* + (a) = vr(±-l) > Oandlm^(x) = tt(±-1) = Im<^_(a). 
Finally, one can show that for x G (— oo, 0), Im 4>+(x) >7r(i — 1) > 0, Im^(x) is decreasing and Re cj)* + (x) 
is increasing. 



Coupling (3.7) and (3.9), we have 



4>*(z)±Tri(l- -) 
c 



from which it is easy to verify the image of the upper half-plane under the mapping eft shown in Figure 
2c. 

To verify the graph in Figure 2d, we first show that for x G (— oo,0), Re^4-(x) <Re</>+(0) < 
and Im.(j) + (x) = —ir. Then, we show that for x G (0,1), Im(f> + (x) < and Im^>_|_(x) is an increasing 
function. To deal with Rec/> + (x) in (0,1), we consider three intervals (0, a), (a, 6) and (ft, 1), separately. 
For x G (b, 1), we show that (Re0 + (x))' = (</>*(x))' < 0. Thus, Re<f> + (x) is a decreasing function. For 
x G (a, 6), we show that Re<^+(x) = 0. For x G (0, a), we first show Re^ + (x) = ^(x). By (3.11) and 
(3.19), Re0+(x) < and (Re<^>+(x))' = {(j){x))' > 0. Thus, Re </>_)_ (x) is an increasing function. This ends 



the verification of Figure 2d, and completes the proof of the proposition. 
To construct our global parametrix, we first define the function 

/ \fz — a + y 1 z — b 

N(z) = 



□ 



2(z-a) 1 /4(z_ 6)1/4 
. v ' z — a — V ' z — b 



V 'z — a — V 1 z — b \ 
2(z-a)V 4 (z- 6)1/4 

y/z — a + V ' z — b 



(3.24) 



V 2(z-o) 1 /4( z _ 6)V4 2 (z - a)V 4 (z - 6)V 4 / 
It is readily verified that N(z) is analytic in C \ [a, b] and 



N + (x) = AL(x) 



1 



x G (a, 6). 



(3.25) 



Next, we introduce the Airy parametrix 



A(z) 



Ai(z) w 2 Ai(u; 2 z) 
zAi'(z) iwAi'(o; 2 z) 



zeC 



+ ■ 



/ Ai(z) 



— ujAi(ujz) 
-ioj 2 A\'(ujz) 



(3.26) 



z G 



where w = e 2?n / 3 . In view of the well-known identity [10, (9.2.12)] 

Ai(z) + wAi(cjz) + w 2 Ai(w 2 z) = 0, 



(3.27) 



10 



it is clear that 



By (3.27), we also have 



A+(z) = A-(z) 



1 -1 
1 



x E 



(3.28) 



k(z) 



1 
±1 1 



—u)Ai(u>z) uj 2 ki(oj 2 z) 
-iuj 2 k\'(ujz) iujk\'({jj 2 z) 



z E 



(3.29) 



-oj 2 Ai(w 2 z) -uki(ojz) 
,—iujk\'{uj 2 z) —iuj 2 k\'(u}z) 

Recall the asymptotic expansions of the Airy function and its derivative [10, 9.7(h)] 



Ai(z) 



-1/4 

1 .2.3/2 



-e s 



E 



(-I)X 



, Ai'(z) 



1/4 00 

' 7 .2-3/2 



E 



(-i) s ^ 



(3.30) 



as z — > oo in |arg z| < 7r, where u s ,v s are constants with no = vo = 1. From (3.26) and (3.30), we obtain 

z - CT3 /4 



A(*) 



2^ 



1 -i 
-i 1 



(I + 0(\z\-'l))e- 2 3 z3/2a:i 



as z — > oo in |arg z| < it, 

k(z) 



1 
±1 1 



-0-3/4 



2^ 



1 



-i 1 



^ (/ + 0(|zrt))e-i 23/2<T3 



as z — > 00 with argz £ (±7r/3,±7r]. Finally, we introduce the D-functions: 

6^(^ + 1/2) 



and 



:= 



2tt(A^) a ^ 
Te Nz ( _ Nz yNz 



(3.31) 



(3.32) 



(3.33) 



(3.34) 



T(-Nz + 1/2) ' 

which were first used in [15] to construct global asymptotic formulas without any cut in the complex 
plane. For simplicity, we also introduce the notations 

D*(z) := D(l -z), D*(z) := 5(1 - z). (3.35) 

By Euler's reflection formula, we have 

D{z) = 2cos(Nirz)e ±N7Tiz D(z) (3.36) 

for z £ C±. The reason why we consider D{z) here is to make sure that the jump matrix Js(z) in Lemma 
4.1 is asymptotically equal to the identity matrix. As n — > 00, by applying Stirling's formula to (3.33) 
and (3.34), we obtain 

D(z) = l + 0(l/n) (3.37) 

for |argz| < it, and 

D(z) = l + 0(l/n) (3.38) 

for |arg(— z)\ < it. 
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i 




TT 

11 









a 


x 




b 


1 


X 




III 






IV 





















Figure 3: The domains I, II, III and IV. 



4 Construction of Parametrix 



First, we select an arbitrary fixed point xo £ Let T be the line Kez = xq; see Figure 3. To 

facilitate our discussion below, we divide the complex plane into four regions by using T and the real line. 
With a similar technique as used in [3] and [15], we construct the parametrix of the RHP for R in these 
four regions. Define 



R{z) := {-l) n ^{We nl y i/2 N{z) 



1 i 
i 1 



/(z)-^/V 1 A(/)a 1 



4 cos 2 (Nkz)D(zY 



N-l 
3=0 



W \[{z- x Nj ) 2 



0-3/2 



for z £ I and III, 

R(z) ■= (-l) N ^(We nl f i/2 N{z) 



\ 1 i V*(z)~^/4 (T2A(/ * )a -l 



4cos 2 (iV7rz)L»*(z) 



N-l 

w n (z-x NJ r 

j=0 



0-3/2 



for z £ II and IV, where f(z) and f*(z) are given by 

2/3 



/(*) 



■-ruj>{z\ 



■-n(j)*(z) 



2/3 



and o"i : = 



1 

1 



,0-2 := 



1 
-1 



(4.1) 



(4.2) 



(4.3) 



. To show that R(z) has the same behavior as R(z) as z — > 00, we 



only consider the case when z S I. The other cases can be handled in a similar manner. A combination 



12 



of (3.31) and (4.3) gives 



-0-3 



(4.4) 



(4.5) 



V^F r. J J /(z)- <J3/4 CJiA(/) C Ji ~ (_l)n e -n(l/2-g)a Se Nniz* s _ 

Here, we have made use of (3.5) and (3.9). Thus, we obtain from (3.36) and (4.1) 

JV-l 

R{z) ~ {We nl Y' il2 N{z)e- n{l/2 ~ g)a:A D^W^ 2 ]J(z- x NJ ) 

3=0 

Also note that D(z) ~ 1 and IT^o 1 ^ ~~ x N,j) ~ z N a s z — > 00. From (3.4) and the asymptotic behavior 
of N(z), it follows that 

R(z) ~ (^e n/ ) CT3/2 Af(z)(^e ni )- <T3/2 e n9CT3 z- Ar ' T3 = (J + 0(l/z))z^ n - N ^ (4.6) 

as z 00. Define the matrix 

S{z) := {We nl )- a ^ 2 R{z)R{z)- l {We nl )' J ^ 2 . (4.7) 

It is easy to see that 

S(z)=I + 0(l/z) (4.8) 
as z — > 00. The jump matrix of S(z) is given by 

J s (z) := S^zy'S+iz) = {We nl )-^l 2 R^z)J R (z)R + {z)-\We nl y*l 2 , (4.9) 
and the contour associated with the matrix Js(z) is S5 = S + T + M; see Figures 1 and 3. 
Lemma 4.1. Js(z) = J + O(^) and S'(z) = I + O (i) as n -> 00. 

Proof. For z € T and ±Imz 6 (0,5), we have Jr(z) = I. This together with (4.1), (4.2) and (4.9) gives 



J s (z) = (-l) N N(z) 
On account of (4.3), we have from (3.32) 



\ *)f{zT^o 2 A{f)o- 2 

-1 



1 i 

i 



(D*(z)/D(z)Y 



Nizrn-ir 



(4.10) 



1 -i 



( ^ / ) /* W^WC/*)^ 1 = (J + 0(l/n))e" 



•n<p 0-3 



1 ±1 
1 



and 



{) /(*)~ a,/4 <^(7)*i 



1 ^ 1 j e ^3 (/ + 0(1/n)) 
as n — )• 00. Applying the above two equations to (4.10) gives 

£>(*) 



J s (z) = JV(z) 



A^(z)- 1 (/ + 0(l/?i)) 



V 

J + 0(l/n) 



D*{z) 
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Here, we have used the fact that D(z) = 1 + 0(l/n), D*(z) = 1 + 0(l/n) and Re^> > 0; see Figure 2c. 
(A corresponding drawing can be given for the image of the lower half-plane under the mapping (p.) 
For z G r and ±Imz ^ (0,(5), we have from (3.31) and (4.3) 

\ ~A r(z)-^ 4 a 2 A(f*)a^ = e-"^(/ + 0(l/n)) 



and 



V5F ( ■ I) Rzy^^AU)^ 



Applying the last two equations to (4.10) yields 
J s {z) = (-l) Af - n AT(z)e- n ^ <T3 
= / + 0(l/n). 



D(z) 



M"(I + 0(l/n)). 



M°*N(z)-\l + 0(l/n)) 



Let us now consider z in the left half-plane of T. For z = x G (— oo,0), we have Jr(x) = I. Coupling 
(4.1) and (4.9) gives 



J s (.) = (-l) n N(z) 



1 i 



1 i 
i 1 



/_(2;) -CT3/4 (tiA(/_)cti 
i 



M uizy^AiU)*! 



{D_(z)/D + {z)r 
N{z)-\-\) n . 



(4.11) 



Note that arg/±(x) G (— 7t,7t) in this case; see Figure 2c. Hence, we obtain from (3.31) and (4.3) 



1 i 

i 



1 



[J f±(z)-^a 1 A(f ± )a 1 = — e - n ^(I + 0(l/n)). 

Substituting the above equation in (4.11) yields 

J s (z) = N{z)e- n ^- a3 e 2Nniza3 e n ^+ CT3 N(z)- 1 {I + 0(l/n)) 
= I+0(l/n). 

Here, we have used the facts D-/D+ = e 2Nniz and e «(?+-?-) = e -Mm*^ 
For z = x G [0, a^o], we have from (2.10) and (2.11) 

/ 1 0\ 



4cos 2 (A^7rx) 
iV-i 

w n (* - ^j) 2 

\ i=o 



/ 



Note that /(z) G C T when z G C±. This together with (4.1) and (4.9) gives 



Js(x) 



N(x)[] I) /(xr^AiU)^ 



1 

D~ 2 1 
-l 



(4.12) 
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From (3.28), it follows that 



Js(x) 



1 i 



1 -1 
1 



0"! 



1 
D- 2 1 



N(x)[] [) /W^V(/ + ) 



1 1 



(4.13) 



Note that N(z) and /(z) x / 4 are both discontinues on the interval (a, xq]. From Figure 2c and (4.3), 
we observe that &rg(f)±(x) = =F7r/2 and arg/±(a;) = =F7r for x G (a, £o]- Thus, we have f + (x)~ as ' 4: = 
/_(x)" CT3 / 4 e i7rCT3 / 2 . By (3.25), it is readily seen that 



N(z) 



1 i 
i 1 



/(*) 



-o- 3 /4 



has no jump on the interval (a,xo\. Applying (3.37) to (4.13) gives Js(x) = I + 0(l/n). 
For z G E-t and Kez < Xq, the jump matrix Jr(z) is 



see (2.12). This together with (4.1) gives 



1 0\ 

- F2e ±i 7 r7Vz cog^TTz) 

^ n (^-^o 2 

V 0=0 1 



J s (z) = N(z) 



1 i 



=Fe 



1 

±iirNz 



\ 



1 i 

i 



\2cos(Nttz)D 2 / 
N-\z). 



[J f{z)-^a 1 A{f)a l 
From (3.31) and(4.3), we then obtain 

j) fizr^Aif)^ = + 0(l/n)). 



1 t 

i 



(4.14) 
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Applying the above formula to (4.14) yields 



/ 



J s (z) = N(z) 



N(z) 



\ 



- s - e 2n(j> e ±2i-KNz 



N{zY l {I + 0{l/n)) 



\2cos(Nitz)e±™ Nz D(z) 2 J 
/I \ 



=Fe 



2n<t> 



N{zY 1 {I + 0{l/n)). 



\D(z)D(z) ) 



One can show that Re</> < in this case; see Figure 2d. Thus, we have Js{z) = I + 0(l/n) 
In a similar manner, we can prove that for z in the right half-plane of T, the corresponding 
Js(z) on £ and on the real line tends to the identity matrix. Hence, we have shown that Js(z) 
as n — > oo on the contour Eg. By the main result in [13], we conclude that S(z) = I + 0(- 



as n — > oo. 
jump matrix 
= /+0(l/n) 
) as n — > oo. 

□ 



5 Main Results 

Theorem 5.1. Let I, D(z) and f(z) be defined as in (3.5), (3.33) and (4-3), respectively. We have 
TT N , n (Nz - 1/2) = {-N) n ^e nl ' 2 

x J r s in(A%z)Ai(/(z)) + cos(Nnz)D(z)Bi(f(z))] A(z,n) 



+ 



sm(NTrz)Ai'(f(z)) + cos{N7Tz)D(z)Bi'(f(z)) B(z, n) 



(5.1) 



for z € I and III, where 



A (*> n ) = {l_a)^ f{z)1/i [1 + ° (1/n)] ' B(Z ' U) = (z- 6)1/4 /(i8) " 1/4[1 + °(Vn)]- 
Similarly, with D*{z) and f*{z) defined in (3.35) and (4-3), 
n N>n (Nz - 1/2) = {-l) N N n ^e nl ' 2 

cos{N7rz)D*(z)Bi(f*(z)) - sin(iV7rz)Ai(/*(z))l A*(«,n) 



+ 



cos(AT7rz)D*(z)Bi / (/*(z)) - sin(JV7rz)Ai'(/*(z)) 



(5.2) 



/or z £ II and IV, where 
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Proof. From the definition of S(z) in (4.7), we have 

R(z) = {We n Y Al2 S{z){We nl )-^' 2 R{z). 
For any matrix X, we denote its element by Xij. The above formula then gives 

Ru(z) = S u (z)R u (z) + S 12 (z)R 21 (z)We nl 

and 

Ri2(z) = S n {z)R l2 {z) + S l2 (z)R 22 (z)We nl . 

First, let us restrict z to the region I indicated in Figure 3. A combination of (3.24) , (3.26) and (4.1) 
gives 



R u ( z ) = (-l) n ^ e nl/2 \-iuM(ujf)a(z,n) - iu 2 Ai'(ujf)b(z,n) 



2cos(Nirz)D{z) 
Ufjo 1 ^ ~ x N,j) 



and 



where 



R 12 (z) = (-l) n V^e nl/2 iM(f)a(z,n) + %M'(f)b(z,n) 



2cos(Nttz)D{z) 



a(z,n) = I" b ^ f(z) 1/ \S 11 (z)-iS 12 (z)) 



and 



(z - a)V4 

= i ^^ f^r 1/ HSn(z)+iS 12 (z)). 

From (2.5) and (2.6), we know that Hn(z) has different expressions in different parts of region I. Let 
us first consider the regions In f2+ and Illn For z 6 I n Q+, we have from (2.5) 

Tie ±N ™ cosjNnz) 
H u {z) = Rn(z) - N _ v R M 

N7T \[j=o {z-x Nd y 

N-l 

= (-l)"^e nl / 2 H(z-x Njj r 1 

x <^ cos(iV7rz)Z)(z) -2iuAi(u)f)a(z,n) - 2ioo 2 Ai'(uf)b(z, n 



- i^ Nz D{z)- 1 Ai(/) a(z, n) + Ai'(/) b(z, n) 
The terms in curly brackets in (5.3) can be rewritten as 



cos(Nnz)D(z) 



■2iuM(uf) -iM(J) -i(D(z)- 2 - l)Ai(/)) a(z,n) 



+ (-2i W 2 Ai , (u;7)-iAi , (7)-i(D(^)- 2 - l)Ai'(/)) b(z,n) 



+ sm(Nirz)D(z)- 1 



Ai(7) a(z, n) + Ai'(7) K z > n ) 



(5.3) 



(5.4) 
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Recall the well-known formula of the Airy functions [10, (9.2.11)] 



Bi(z) = ±i 



2e^/ 3 Ai(; ^ 



CO z) 



A\{z) 



(5.5) 



Also, note that D(z) ~ 1 for z / and Ai(/) and Ai'(/) are exponentially small as n — > oo when z is 
close to the origin. Hence, we can always neglect the terms (D(z)~ 2 — l)Ai(/) and (D(z)~ 2 — l)Ai'(/) in 
(5.4). This together with Lemma 4.1 and (5.5) gives 

N-l 

H n (z) = (-l) n ^e nl ' 2 H (z - x^-)" 1 

3=0 

sm(NTTz)Ai(f(z)) + cos(NTTz)D(z)Bi(f(z))] A(z,n) 



+ 



sm(N7rz)Ai'(f(z)) + cos(N7rz)D(z)Bi'(f(z)) B(z, n) 



Similarly, for z G III H J7_ 



JV-i 



j=0 



x <( cos(Nttz)D(z) 2iu 2 Ai(oo 2 f)a(z,n) + 2iuAi' (u z f) b(z,n 



i i. .2 : 



(5.6) 



(5.7) 



+ ie-^ Nz D{z)~ 1 Ui(/) a(z, n) + Ai'(/) b(z, n) 



(5i 



Again, by (5.5), we get exactly the same formula given in (5.6). Using Lemma 4.1, we can obtain the 
desired result (5.1). Now, we show that the asymptotic formula of Hn(z) in (5.6) holds not only for z 
in £l±, but also for z in the whole I and III. From the relation between R(z) and H(z) in (2.6), we know 
that H\i{z) =R\i(z) for z G I In contrast to the expansion in (5.3), for 2 6 I \f2+ the term 

N-l 

{-l) n ^e nl ' 2 H(z- xn^)- 1 x { - ie^'Diz)' 1 [M(J) a(z, n) + Ai'(/) b(z,n)] } 
3=0 

does not appear, since the quantity e lnNz Ai(f) is exponentially small as n goes to infinity, in comparison 
with the other term in (5.3). This suggests that the region of validity of the expansion given in (5.6) can 
be extended to z G I U III. As a consequence, (5.1) holds for z G I U III. 
In a similar manner, we have 

N-l 

H u (z) = (-l) N V^e nl / 2 ]](z- x^-)" 1 

3=0 

x J cos(Nirz)D*(z) [2ioj 2 Ai(uj 2 f*)a*(z,n) + 2iujAi(u} 2 f*)b*(z,n)] 



+ ie™ Nz D*(zY l [Ai(/*) a*(z, n) + Ai'(/*) b*(z, n)] 
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where 

= 7^^ /*(^) 1/4 (gll+^12) 

and 

b *^ n ) = { "~ b) 'l r^r l, \s ll -is 12 ) 

(z — a) 1 / 4 

for z £ II H f2+. Following the same argument as given above, one can obtain (5.2) for z E II and IV. This 
completes the proof of Theorem 5.1. □ 



6 Comparison with Earlier Results 

To compare our formulas in Theorem 5.1 with those given in [12], we first derive from (5.1) two simple 
asymptotic formulas for TT^^ n {Nz — 1/2) when z is real and less than a; cf. (3.2). 

Theorem 6.1. Let I, 4>{z), D{z) and D(z) be defined as in (3.5), (3.9), (3.33) and (3.34), respectively. 
We have 

ir N , n (Nz - 1/2) = {-NYe nl ' 2 {e-^D{z) cos(iV^) ^-^+i^Z^! [1 + (l/n)] 



+ 0(e nRe ^ z) )j (6.1) 

for < z < a — 5 < a, and 

n N , n (Nz - 1/2) = Nnenl/2 {^ nHz) D(z) { l~ a ^ [1 + 0(l/n)} + 0(e nRe ^)} (6.2) 

for z < 0. 

Proof. From the well-known asymptotic expansions of the Airy function [10, (9.7.5)-(9.7.8)], we have 



Ai(7<») ~ f; -t^r, Ai'(7 ( ,)) ~ JMIt^M f; m 

2 ^ 2 ^ ^(-n^>(z)) fc 



and 



DC 



Bi (7(,)) ~ ZW_^ e -«*w y — ^ — Bi'(7(.)) ~ M^ e -«^) y — v j — (6 . 4) 

where /(z) is defined in (4.3). Note that <p{z) is analytic for real z 6 (0, a); see the paragraph containing 
(3.9). Applying (6.3)-(6.4) to (5.1) gives 

n N , n (Nz - 1/2) = {-N) n e nl / 2 j e »?M sin(jV^) ^^^"'^ t 1 + 0(V*)] 

+ e-**>2>(*) cos(iVvr.) ^f^^t [1 + 0(V«)] } ■ (6-5) 
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Since (j>(z) < in this case on account of (3.11), e n ^ z ^ is exponentially small as n goes to infinity. Thus, 
(6.1) follows from (6.5). 

Note that 4>(z) has a jump across the negative real axis. For real z < 0, we obtain from (6.3)-(6.4) 
and (5.1) 

7T+ n (Nz - 1/2) = (—N) n e nl / 2 jc^W MN^z) ^^^^^ [1 + 0(l/n)] 

+ e-^Wfl + (z) cob(Mtz) ( ";_t)V^ir &) V4 /2 [1 + °(Vn)]} (6.6) 

as n — > oo, where ir^ (Nz — 1/2) denotes the limiting value of iiN,n{Nz — 1/2) as z approaches the real 
line from above. Observe from Figure 2c that Re0+(z) < 0. Thus, it follows from (6.6) that 

n+jNz - 1/2) = (-AOV^je-^W D+(z) cos(Nnz) {z ~ ^ ^ ~*f'* [1 + 0(l/n)] 



(z-a)V4( z _ 6)V4 



+ 0(e 



nRe </>(z) ' 



(6.7) 



Similarly, one can see that it n n (Nz — 1/2) is given by 

(-ty» e «/*[ e -M-W D _( z) cos(iVvrz) ^Z^+i^Z^! [i + 0( i /n)] + 0(e nBe? W) J_ (6 . 8) 

From the definition of in (3.5), we note that e n ^ z ^ can be analytically extended to (— oo, 0). Thus, from 
(3.9) and (3.36) it follows that 



e -n*±(.*)D ± (z)coB(Nirz) 



(-ir 



~ n ^D(z). 



(6.9) 



In view of (6.9), the two asymptotic formulas (6.7) and (6.8) are exactly the same. Hence 

\i/2 + ( x -hW 2 



KN,n( Nz ~ V 2 ) = N n e nl/2 
for real z < 0. This gives (6.2). 



(6.10) 



□ 



We will now show that our asymptotic formulas for the discrete Chebyshev polynomials are the same 
as those given by Pan and Wong [12]. First, we introduce the notation 



x := Nz - 1/2. 



(6.11) 



Two different asymptotic approximations for t n (x,N + 1) are given in [12, (8.13) and (8.6)], one for x 
negative and the other for x positive. Changing N + 1 to N, they read 



tn{x,N) 



l) n+1 r(n + N + 1)(N- l) x n- 2x - 2 T{x + 1 



x < sin(7rx) 



T(N)tt 



1 + 0| iv 



+ 0{e 



(6.12) 
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for fixed x > and 



tn(x,N) 



(-l) n r(n + N + 1)(JV - l) x n 
T(N)T(-x) 



x„-2x-2 



1 + 01 N 



(6.13) 



for fixed x < 0. 

We next derive from (6.1) and (6.2) asymptotic formulas for t n (x, N) when x is fixed (i.e., z = 0(1/N)). 
For x > 0, i.e., 2 > 0, substituting (6.11) in (3.33) gives 



e Nz T(Nz + 1/2) 



D(z)cos(Nirz) 



2it(Nz) Nz 



cos(Nnz) 



e x+l/2 T {x + 1) 



2ir{x + l/2) x + 1 l 2 



sin(7rx). 



Moreover, it is readily verified that 



|2 



7T 



as n — )• 00 and 

( Z - a) 1 /* + (* _ &)V2 a l/2 + 6 l//2 ( a + fc + ^1/2^,1/2)1/2 



(2 - a)V4(z - 6)1/4 ( a6 )i/4 



(a6)V 4 



V2(l + c) 1 '\ 



-1/2 



(6.14) 



(6.15) 



(6.16) 



as z — > 0. Here, we have made use of the fact that a + b = 1 and afe = c 2 /4; see (3.2). 
Note that z — > as n — > 00; a combination of (2.1), (6.14)-(6.16) and (6.1) yields 

t n (x, N) ~ ( _i ) ^+i iV »2 2 "n- 1 /2 J^^^Xvl (1 + c) 17 ^" 1 / 2 sin( 7 rx) e "( i / 2 -^» 



(6.17) 



as n — > 00. We now derive an explicit formula for 1/2 — 4>{z); see (6.26) below. From (3.5) and (3.9), one 
has 



1 



1/2 - <j){z) = 1/2 - <f) + (z) - m{\ - -z) = Re g+{z) 



(6.18) 



for < z < a — 5. This fact evokes us to calculate g{z) first. From (3.3) and (3.4), it is easily seen that 



9{z) 



1 ra 1 r 1 

— / log(z — s)ds H — / log 
c Jo c Jb 



By integration by parts twice, we obtain 

rb 



log(z — s)fi(s)ds = fi(s) log(z — s)s — zfi(s) log(z — s) 



(z — s)ds + / log(z — s)/j(s)ds. 



fi(s)ds 



(6.19) 



pb pb 

+ z fi(s)' log(z — s)ds — / s log(z — s)fi(s)'ds. 

J a J a 



By the definition of fi(x), we have 



1 - 2x 



2vr (x _ x 2 )^/(x - a)(b - x) 



(6.20) 



(6.21) 
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Here, we have again used the relations a + b = 1 and ab = c 2 /4. Moreover, one can show that 



log (2; — s)fi(s)ds 



-(z — a) log(z — a) (z — b) log (2 — b) 

c c 

b log(z-s) 1 



+ 



2vr7 a Sv /( s -a)(6- 5 ) 
log(z - s) 



2vr 



log(z - s) 



a (l-s)V(«-a)(6-s) 



ds 



ds. 



(6.22) 



Let the integrals on the right-hand side of the equality be denoted by I±, I2 and ^3, respectively. To 
evaluate I±, we note that from [16] we have 



log s 



zds 



2tt 



L (s-zW(s-a)(b-s) (z-ay/2(z - 6)V2 
Making the change of variable s = z — x, (6.23) gives 

logx 



log 



z-b 



+ V^b+(z-a) 1 / 2 (z-b) 1 / 2 
(^+Vb)z 



d(z — x) 



(6.23) 



1 (z — x) y/ (z — x — a)(b — z + x) 

4vr , z + J{z-a)(z-b) + c/2 

= log , . 

c (y/z - a + \J z - b)z 

Similarly, we can evaluate I2 and I3. Thus, a straightforward calculation shows that 

g( z ) = -1 - 2 log 2 + -(z - l)log(*- 1) - -zlogz+ (2 - -)log \(z - a) 1/2 + (z - 6) 1/2 
c c c L 



(6.24) 



2 , 

H — zlog 
c 



z + (z - a) l l 2 {z - b) 1 ' 2 + c/2 + -(1 - z) log 2 - 1 + (z - a) 1 / 2 ^ - b) 1 ' 2 - c/2 

J c L 



(6.25) 



By (6.18) and (6.25), 1/2 — (ft(z) can be explicitly given by 

- 1 - 2 log 2 + -(z - 1) log(l - z) - -2 log z + (2 - -) log \(a - z) l/2 + (b - z) 1/2 
c c c L 

2 



H — z log 
c 



z-(a-z) 1 / 2 (6-z) 1 / 2 + c/2 + -(l-z)log I- z + {a- z) 1 ' 2 ^- zfl 2 + c/2 . (6.26) 



Note that n = cN and z — > as n — > 00. From (6.26), we have 



1/2 - cj>{z) 



-l-21og2 + (l + -)log(l + c) 

c 



+ z 



- log(z) log c 

c c c 



+ 0(z 2 



and so 



D n(Z/2-#*)) 



e ~n 2 -2n (1 + c )n+JV (x + ^^l^g-s-l^z+^-to-l^ 



Substituting (6.28) in (6.17) yields 



t n (x,N) 



(-i) n+1 r(x + i) 



-2x-2 



■ sin 



7T 



(^x)iV n+X+1 (l + c )"+^+l/2 e -n 



(6.27) 



(6.28) 



(6.29) 
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which will agree with (6.12) to leading order if we can show that 



T(n + N + l)(N-ir Nn+X+1/ 
T(N) 



n+N+l/2 e -n 



(6.30) 



The last result is a direct consequence of Stirling's formula. 

Next, we consider the case when x is negative, but for a moment let's restrict — oo < x < — \ (i.e., 
z < 0). Inserting (6.11) in (3.34), we have 



D(z) 



2^e x+l ' 2 



(-x-l/2) x + l / 2 r(-x)' 
A combination of (6.2), (6.15)-(6.16) and (6.31) gives 

e x+l/2 



t n (x,N) ~ N n 2 2n n~^ 2 



(-x- l/2)*+ l / 2 T(-x) 



n +c )l/2 c -l/2 e n(V2-<!!.(2))_ 



(6.31) 



(6.32) 



Since e n< ^ 2 ) is analytic in the interval (— oo,0), we obtain from (3.5) that e n ^l 2 = e n ^l 2 
e n a+i z ) _ On account of (6.25), it can be shown that 



9+[z) 



iri - 1 - 2 log 2 + (1 + -) log(l + c) 

c 



+ 2 



- log(— z) log c 

c c c 



+ 0(z 2 



as z — ?• 0. Thus, letting n — > oo, we have from (6.33) 



n(l/2-4>(z)) 



4) n 2" 2n (l + c) n+7V e- n C - 2a; - 1 ^^ x ^ 1/2 (-a; - l/2) x+1 / 2 e - x - 1 



/2 



From (6.32) and (6.34), we obtain 



tn(x,N)~(-N) n 



n 



-2x-2 



V{-x) 



( l + c) n+N+l/2 N x+l e -n_ 



(6.33) 



(6.34) 



(6.35) 



In the case — i < x < 0, formula (6.35) follows directly from (6.29) on an appeal to the reflection formula 
r(l + x)T{— x) = — it/ sin(7rx). On the other hand, coupling (6.13) and (6.30) gives 



t n (x,N)~ (-JV)' 



, n 



-2x-2 



T(-x) 



(1 + c) n+N+l/2 N x+l e ~n 



(6.36) 



as n — > oo. Therefore, our results (6.29) and (6.35) agree with those of Pan and Wong [12] stated in (6.12) 
and (6.13). 
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